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Abstract 

In a recent paper, Jones and Moore have deve- 
loped a simple numerical lifting surface technique 
for calculating the aerodynamic coefficients on 
oscillating wings in subsonic flight. The method 
is based on the use of the full lifting surface 
theory and is not restricted in any way as to fre- 
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M < 1. In this study, this simple but general met- 
hod of predicting airloads is applied to helicopter 
rotor blades on a full three-dimensional basis. 

/The general theory is developed for a rotor blade 
at the ip « tt/ 2 position where flutter is most lik- 
ely to occur. Calculations of aerodynamic coeffi- 
cients for use in flutter analysis are made for 
forward and hovering flight with low inflow for 
Mach numbers 0 and 0.8 and frequency ratios p/ft* 1 * 
and 4. The results are compared with values given 
by two-dimensional strip theory for a rigid rotor 
hinged at its root. The comparisons indicate the 
inadequacies of strip theory for airload predicti- 
on. One important conclusion drawn from this stu- 
dy is that the curved wake has a substantial effect 
on the chordwise load distribution. The pitching 
moment aerodynamic coefficients differ appreciably 
from the results given by strip theory. 

Introduction 

In a survey paper. Ref. 1, Jones' et al. give 
a detailed account of significant developments in 
the field of unsteady aerodynamics of helicopter 
rotor blades. One of the problem areas surveyed 
was that of blade flutter as it has been found that 
under certain operating conditions, rotor blades 
can flutter in both hovering and forward flight. 
This phenomenon has been investigated by several 
researchers in Refs. 2, 3, 4, and 5 and the results 
of their studies have improved our understanding 
of the problem. For the case of hovering flight, 

J. P. Jones in Ref. 2 applied a method developed 
by W. P. Jones in Ref. 6 to derive the approximate 
aerodynamic coefficients for an oscillating single 
rotor blade for use in his flutter analysis. He 
approximated the actual flow conditions by neglect- 
ing curvature effects and assuming a simple two- 
dimensional mathematical model cosisting of a ref- 
erence blade and an infinite number of wakes lying 
beneath the reference blade extending from - 00 to ®. 
He considered flapping and pitching motions and com- 
pared his results with those obtained experimenta- 
lly by Daughaday and Kline in Ref. 3. On the basis 
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of this work it was concluded that the wake is pri- 
marily responsible for some of the vibratory pheno- 
mena found on helicopters in practice. For low 
inflow conditions, Loewy in Ref. 4 used a similar 
mathematical model to that of J. P. Jones and inve- 
stigated the variation in the pitching moment damp- 
ing coefficient of a particular blade section as 
p/fl varied for specified positions of axis of osci- 
llation and a range of values of wake spacing. He 
found that the damping coefficient became negative 
whenever p/ft was slightly greater than an integer 
for axis of oscillation forward of quarter-chord. 
Similarly he found that the d amp ing coefficient 
for a flapping oscillation dropped sharply at inte- 
gral values of p/ft but did not actually become neg- 
ative. Timman and Van de Vooren in Ref. 5, on the 
other hand, assumed that there was no Inflow thro- 
ugh the rotor disk and developed a theory for cal- 
culating the aerodynamic forces on a blade rotating 
through its own wake. Their results agree with 
those obtained in Refs. 2 and 4 in the limit when 
zero spacing between the wakes is assumed. All 
this theoretical work confirms the conclusion that 
the proximity of the wake is a contributing factor 
to rotor blade flutter. 

All the theoretical work described above is 
based on the assumption that the flow is incompre- 
ssible. However, with the advent of helicopters 
capable of flying with blade tip speeds ranging up 
to and in excess of the speed of sound, compressi- 
bility effects need to be taken Into account when 
determining coefficients for use in flutter analy- 
sis. Jones and Rao in Ref. 7 were able to do this 
on the basis of two-dimensional theory and have 
derived coefficients for a range of Mach numbers, 
reduced frequencies, and wake spacing. Their ana- 
lysis is based on the use of Loewy's model. Ref. 4, 
of the helical wake and the application of a theory 
developed earlier by Jones in Ref. 8 for an oscilla- 
ting airfoil in compressible flow. The values of 
the coefficients given in Ref. 7 agree with those 
obtained in Refs. 2 and 4 for zero Mach number but 
differ appreciably when the Mach number is varied. 
Hamnond in Ref. 9 also developed a theory for det- 
ermining compressibilty effects by using a differe- 
nt model of flow from that used in Ref. 7. In his 
model, the wake of the qth blade of a Q bladed 
rotor after n revolutions extends from — 2tt (n+q/Q) 
to ®; in Jones and Rao's model it extends from -°° 
to His aerodynamic coefficients for several 
Mach n umb ers and inflow ratios are in general agr- 
eement with* the results of Jones and Rao in Ref. 7. 

While the aerodynamic derivatives predicted 
by two-dimensional strip theory are widely used in 
predicting the flutter speeds of helicopter rotor 
blades , the method does not allow for curvature 
and finite aspect ratio effects. For incompressi- 
ble flow, Ashley, Moser, and Dugundji in Ref. 10 
developed a three-dimensional model in which they 
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modified Reissner's theory. Ref. 11, for oscillat- 
ing wings in rectilinear flow by including the free 
stream velocity variations along the span. Their 
results indicate a negligible difference between 
two and three-dimensional solutions up to 95% of 
the span. Jones and Rao in Ref. 12 similarly stu- 
died tip vortex effects in compressible flow and 
they also concluded that such effects are negligi- 
ble except in regions close to the tip. In some 
of his earlier work. Miller in Refs. 13, 14, and 
15, developed a helical wake model in which the 
rotor wake was divided into a "near" wake and a 
"far" wake. The near wake included the portion 
attached to the blade that extend approximately 
one-quarter of a revolution from the blade trailing 
edge. The effects of the near wake include an in- 
duced chordwise variation in downwash and were for- 
mulated using an adaptation of Loewy's strip theory. 
The chordwise variation in the velocity over the 
airfoil induced by the far wake was neglected. 

Miller extended his model to study the forward fli- 
ght case and found that the nonuniform downwash 
induced at the rotor disk by the wake vortex system 
could account for the higher harmonic airloads en- 
countered on rotor blades in forward flight. He 
also showed that under certain conditions of low 
inflow and low speed transition flight the return- 
ing wake could be sucked up into the leading edge 
of the rotor which would account for some of the 
vibration and noise. Piziali in Ref. 16 has deve- 
loped ^n alternative numerical method in which the 
wake of a rotor blade is represented by discrete 
straight line shed and trailing vortex elements. 

He satisfied the chordwise boundary conditions, but 
the rotor blade was limited to one degree of free- 
dom in flapping. Sadler in Ref. 17, using a model 
similar to Piziali* s, developed a method for predi- 
cting the helicopter wake geometry at a "start up" 
configuration. He represented the wake by a fine 
mesh of transverse and trailing vortices starting 
with the first movement of the rotor blade genera- 
ting a bound vortex, and, to preserve zero total 
vorticity, a corresponding shed vortex in the wake. 
Integrating the mutual interference of the trailing 
and shed vortices upon each other over small inter- 
vals of time, Sadler was able to predict a wake 
geometry. Although his model showed fair agreement 
with the available experimental data for advance 
ratios above one-tenth, Sadler’s method is limited 
due to the large computational time required tp 
represent the wake by a finite mesh. 


soning led them to represent the blade motion by a 
series of oscillatory pulses, where each disturba- 
nce occurs over the range, y - < i|; < y + 

Corresponding to each burst of oscillation, packets 
of vorticity are assumed to be shed into the wake. 
With increasing forward speed, the spacing between 
the packets of vorticity also increases and it was 
found that the flutter speed became constant when 
y, the advance ratio, was above 0.2. The approach 
used in the present study differs from that adopted 
by Shipman and Wood in that continuous high freque- 
ncy small oscillations are assumed to be superimpo- 
sed on the normal periodic motion of the blade. 

The airloads and aerodynamic derivatives associated 
with the perturbed oscillation of the rotor blade 
can then be calculated by the method described in 
this paper. Since the rotor blade will first att- 
ain its critical speed for classical flutter at 
® tt/2, the aerodynamic derivatives corresponding 
to this value of iJj only have been calculated. The 
method takes finite aspect ratio and subsonic com- 
pressibility effects fully into account. Typical 
results for a rotor blade hinged at its root desc- 
ribing flapping and twisting oscillations are given 
for a range of Mach numbers and frequency values. 

Basic Equations 


In the development of the analysis of the 
Jones-Moore theory, Ref. 18, for oscillating wings 
in rectilinear flight, the space variables x, y, z, 
and t are replaced by X, Y, Z, and T, respectively, 

so that- IZ . IT m 

x £X, y “ ^ > z “ g > t - u 

where £ is a convenient reference length, U is 
the uniform velocity, M is the Mach number and 
2 1/2 

8 *= (1-M ) 7 . The velocity potential of the flow 
around a surface oscillating at a frequency p 
can then be expressed as 


<Kx,y,z,t) = U£$(X,Y,Z)e 

2 2 
m w/e 


i(XX + wT) 


( 2 ) 


where u) = pJl/U, X = M w/8 . The function $ may be 
regarded as a modified velocity potential. Fur- 
thermore, it can be shown that it satisfies the 
wave equation 
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where k = Mo)/ 8 


Though many forms of flutter can occur on ro- 
tor blades, attention in this report is concentra- 
ted on the determination of appropriate aerodynamic 
coefficients for use in the analysis of blade flu- 
tter of the classical bending-torsion type. Ship- 
man and Wood in Ref. 20 have considered this prob- 
lem but they did not take compressibilty and fini- 
te span effects into account. The two-dimensional 
mathematical model used is similar to that employ- 
ed by other authors except that they assumed that 
flutter would first occur when the relative velo- 
city over the rotor blade reaches its critical va- 
lue when ^ - tt/2. For greater or lower values of 
ij/, the relative speed would be reduced below the 
critical speed for flutter and any incipient grow- 
ing flutter oscillation would be damped. This rea- 


Since in this problem the motion of the sur- 
face is assumed to be prescribed, the downwash 
velocity at any point on it must be the same as 
the downwash induced by the velocity (or doublet) 
distribution over the surface and its wake. This 
condition must be satisfied in order to ensure 
tangential flow over the surface at all points. 

It Is also assumed that the rotor blade is a thin 
surface oscillating about its equilibrium position 

in the plane z = 0. If £ = £ f e* pt defines the 

downwash displacement at any point (x,y) at time 
t, this boundary condition requires that the down- 
ward velocity and 3<J>/3z must be equal. In the 
transformed coordinates, this implies that 


128 


I 


w « 


3* 


-i(XX + i*>T) 
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where w * is known, 

ot ox 

A further condition that must be imposed on 
any solution is that it leads to zero pressure dif- 
ference across the wake created by the oscillating 
surface. From the general equations of flow it can 
be established that the local lift £(x,y,t) at any 
point is given by 

I(x,y,t) - p(|£ + 0 |“) (5) 

where k * <j> u - the discontinuity in the veloci- 
ty potential. From Eq. (5), it immediately follows 
that on the surface 


£(x,y,t) 

2 

where v * co/B 
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and K 
« 0 


This yields 


( 6 ) 
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everywhere in the wake since the lift must then be 
zero. From Eq. (7), it can be deduced that at any 
point in the wake 

K(X,Y) - K(X t ,Y) e _iv(X_X t ) 


( 8 ) 


where X « X^_ denotes the position of the trailing 
edge of the section at Y. 


As shown in Ref. 19, the solution of Eq. (3) 
may then be derived from the integral equation 


4irW(X ,Y ,0) 
P P 
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n K (X , Y) — ) dXdY 
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where W is the modified downwash at the point 

given by Eq. (4), K has to take the form specified 
by Eq. (8) at points in the wake and 

€ = [(X-X p ) 2 + (Y-Y p ) 2 + Z 2 ] 1/2 . 

The double integral in Eq. (9) must be taken 
over the area of the oscillating surface and its 
wake. It should be remembered, however, that K * 0 
along the pleading edge and the sides of the area of 
integration. 


In the numerical technique developed in Ref. 

18 for calculating the airloads on oscillating 
wings in rectilinear flight, the wing is divided 
into a number of conveniently shaped boxes and K 
is assumed to be constant over each box. The wake, 
on the other hand, is divided into a number of cho- 
rdwise strips and K over each strip is defined by 
Eq. (8). The contribution of the wake to the down- 

wash W(X ,Y ) is then derived by direct numerical 
p p 
integration. 

The application of the method outlined above 
to determine the airloads on rotor blades presents 
certain difficulties, the principal one being that 
the flow velocity over the rotor blades is not con- 
stant as assumed in the derivation of Eqs. (3) and 
(8) for wings in straight flight. To overcome this 
difficulty, it is assumed that the rotor blade can 
be represented by a number of spanwise segments 


over every one of which the flow is taken to have 
its average value and appropriate Mach number. On 
this basis the above analysis can be modified for 
application to rotor blades as outlined in the 
next section. 

Rotor Blade Theory 

In the present analysis, the rotor blade is 
taken to be fixed at the $ - tt/2 position and its 
helical wake is assumed to extend rearwards as 
indicated in Fig. 1. Normally, one would expect 
the vorticity shed by the perturbed blade to be 
carried downstream by the distorted wake of the 
loaded rotor blade. However, in the present preli- 
minary study, uniform inflow is assumed and any 
distortion of the wake due to blade-tip vortex 
interference is ignored. The aerodynamic coeffic- 
ients corresponding to any prescribed motion can 
then be calculated for forward and hovering flight 
by the method described below. 

a) Forward Flight (Rotor Blade at $ = ir/2 ) 

Let R denote the tip radius of the blade and 
assume x * Rx* , y ■ Ry* , and z *= Rz f . For forward 
flight with velocity V, the relative local velocity 
at section y will be denoted by U(«V-K2Ry f ) and U 1 

1 ) , where ft is the angular rotation and u(*V/ 
ftR) is the advance ratio. It then follows that at 

the section y f , the downwash w(x f ,y ? ) is given by 
P P P 

w ( x p ,y p> “ QR(i ^’ + d p& eipt (10) 

where z, * Rc*e ipt is the displacement of the blade 
at the point, (x^y^). When the blade is describ- 
ing flapping and twisting motions, £* may be expr- 

peepH OQ 

c' - Y’f(y') + x'a'F(y') (11) 

where y 1 and a* are the amplitudes at the reference 
section and f(y f ) and F(y ! ) are the modes of flap- 
ping and twisting oscillations, respectively. If 
the blade is assumed to be rigid and hinged at the 
root, f(y f ) * y f , and F(y ? ) * 1 in the above equa- 
tion. For convenience, the reference section is 
taken to he at the tip but, in actual flutter cal- 
culations, the section at 0.8R would be a better 
choice. 


To obtain the distribution of R corresponding 
to the motion prescribed by Eq. (10), Eq. (9) is 
first expressed in terms of the original variables 

and K is replaced by Rk , e ipt . It may then be 
written as 


-iX’x* 
4irw f e p p 


//k* (x* ,y’)e 

dx’dy * 
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-iX ' x'3^ ,e~ ilC ' r ' 


where k * 


e 2 nu' 


, X' = Mk' , w - w'e ipt , g 2 


and r’ - [(x'-x^) 2 + g 2 (y’-y^) 2 + g 2 z ,2 ] 1/2 . 


( 12 ) 

1-M 2 , 


The above equation can be used to obtain the 
solution to the problem of determining the flow 
over a rotor blade with a rectilinear wake. Since 
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the wake can withstand no lift, the condition, 
3 k »k 

— + = °, must 

wake, this yields 


3 k 

= o, must be satisfied. For a rectilinear 


-i 


p(x'-x^) 


k' (x',y') - k^(y')e 


fiU' 


(13) 

However, if the wake originating from a blade strip 
is assumed to be curved 


-i 


P(s'-sp 


ftq ? 


k f (s*,y’) = k'(y')e ^ , (14) 

where q f « (y 2 + y ' 2 + 2yy , sin i|*) 1/2 , s « Rs* is 
the distance along the vortex path and y* specifies 
the spanwise location of the blade strip. 

For computational purposes, Eq. (12) may be 
conveniently expressed as 

n. 

(15) 


4irW(x',y') - //K(x'.y') £ dx'dy'. 


“p’'p v ' 8 

_ w(x'y’) -iX'x' 

where W(x^,y^) - 1 — E- e p p , 

P 

K(x',y’) - k'(x',y') e~ iX ' x ’, and 
3 2 e -iK'r' o -iie'r' 

G 2 ( 7^ m <"V ' H < 1+ iic'r') 

32 1 1 

(l-^ 2 li) + k , 2 b 2 z ’ 2 ]. 

r f 2 

It should be noted that in the wake 

V’(8’-S') 

-i 

K(x',y') - K t (y') e q , (16) 

where v' - -E— . 

ne 

b) Hovering Flight 

For the simplest case of hovering flight, 

V - 0 and s' - y'0. Hence Eqa. (10), (11), (12), 
(13), and (15) can be simply modified by replacing 

V with zero. Eqs. (14) and (16) then become 

-if(e-e ) 

k'O.y') « k’(y') e n * (17) 

and - - -iv'(e- 6 .) 

K(9,y') - K'(y') e C , (18) 

where y' defines the location of the blade strip 
from which the wake originates. 

Method Of Solution 

The schematic diagram of the oscillating rotor 
blade is shown in Fig. 1. Eqs. (15) and (16) are 
combined and expressed as 

-ixV 

47 TW(x',y') - - f f K(x\y')B(- r— )(l+iK'r') 

p p blade r ' 3 8 

surface 

X V Z t B t ; -lx’r^ 

dx'dy' - / / K (y’)e ^ B(~ 5 — ■) 

wake r * J 
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2 2 

[ ( 1 +iK’rpd- # z j ■ ) + k’Vz ' 2 ] ds'dn', (19) 

r f 

w 

where r^ = [ (x'-x ') 2 + 8 2 (y'-y ') 2 J 1 ^ 2 

r; = I (x’-xp 2 + B 2 (y'-yp 2 + e 2 z’ 2 ] 1/2 
ds' - (dx' 2 +dy ' 2 ) 1/2 

and dn f is perpendicular to ds* and approximately 
equal to dy f on blade. The rotor blade is divided 
into a number of rectangular boxes (M x N) on which 
the doublet strengths are assumed to be constant as 
in Ref. 18. Based on this assumption, Eq. (19) can 
be expressed as 

= E E S.. K + E T. K (20) 

“ i-i j«i ij lj j-i J tJ 
In Eq. (20) and T^ may be interpreted as the 

aerodynamic influence coefficients and the actual 
expressions are given later in this section. 

Tj are the downwash velocities induced at the box 

ran due to the unit strength doublets located^ at the 
box ij and the wake strip j respectively. and 

are the doublet strengths at the box ij and the 

trailing edge of the wake strip j , respectively. 
With the use of the wake boundary condition, K^_ 
can be expressed as 

1 u 1 


tj 


J 


K tj " 

Eqs. (20) and (21), then yield 

4ttW - E Z A. 4 K 44 
““ i-i j-i 13 13 

where A^ - for i i 1 M 

)'fv' -x' 


v' (x',-x' ,) 

+ 21 ^ *& ■] ( 21 ) 


and A - S ±j + Tj/fe 
for i • M. 




+ 21 - 


( 22 ) 




U i 


For a given mode shape (W f s are known), Eq. (22) 

r p-n 

represents a system of H x N linear algebraic equa- 
tions, the solution of which yields the values for 

K . 's. M and N denote the total number of chord- 

mn 

wise and spanwise stations respectively. 

Once the appropriate K distribution has been 
found, it is then relatively easy to determine the 
aerodynamic forces per unit span acting on the rotor 

blade. If, In Eq. (5), k - Rk'e ±pt - KKe lpt e 1X ' x ' , 

it then follows that the local lift L(y) * L 1 (y) 

e^ t and the nose-up pitching moment, M(y) ■ M f (y) 

e pt , referred to the mid-point of the chordwise 
section at y are given by 


JlIIiL 


P (U R ) (c/ 2 ) 


2 *'t 


(Lj+il^) (-±) + a 3 +lL 4 )o' 
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/ k’dx' + O'k’) 
» *'■ 


(23) 


2 z! 


■7") + <M 3 +iM 4 )a' 


“ (M,+iM,) ( 

P(D_ ) Z (c/2) \ 

*■ t: t 

»(^)(— ) 2 [0’C / k'dx’ - k'x*) - l£ / k'x’dx*] 

9 c 

R i. 4 (24) 

where c Is the local chord, U R Is a reference velo- 
city, z^ and a’ are the local amplitudes of the 
flapping and twisting motions 







ly, and L_ 


L 3 , M^, H 3 , and L 2 , L^, M 2# M^, are the in phase 
and out of phase airload coefficients, respectively. 


Expressions for Aerodynamic Influence Coefficients 

Forward Flight (Rotor Blade at it> = tt/2) 

The influence coefficients are calculated by 
the method outlined in Ref. 18. For a box not con- 
taining the collocation point 


n! 4 d« „ * /« x -iK’r 1 

j 2 « -iv’(e-e ) w 

T. * - / /« z M- r— HU+ifc’r’) 

3 »j-< 2 •« *: 

(1- M V 2 ) + K^fl 2 *' 2 ] 7l d&dn’ (27) 

r ,Z J 

w 

where x* * y’ sin 8 , y* * y! cos 9, and z' « d , 9/2ir. 
j J 

The effect of the helical wake in hovering 
flight is estimated by two different methods. In 
the first method, a Helical Wake Model is used and 
the actual helical path is taken in evaluating the 
Tj coefficients. In the second method, a Circular 

Wake Model is employed and the helical wake is re- 
placed by its near wake, which is assumed to extend 
over © t £ 9 $ ir/ 2 , and a number of regularly spaced 

circular disks of vorticity below the reference 
plane. The formula for k f for the n’th disk at 
z 1 * nd* is taken to be simply 

-if [(9-0 )+2nir] 

k'O.y’.nd') - k’(y')e ( 2 8 > 

the actual spacing between consecutive disks being 
Rd\ 


3 ±i 


^2 

- - I 

y j " d 2 


where ic* * ■ 
d^ « Ax*/2, 


Mp 

6 2 «U* 






-ix'r; 

(r 5 — jea+lK’rpdx'dy’ (25) 

r’ J 8 

s 


, r' - I(x’-x ’) 2 + 6 2 (y’- y ;) 2 ] 1/2 , 
s m n 

Ay'/2, and Ax' and Ay* are the 


chordwise and spanwise s pacings of the rectangular 
grid on the surface of the rotor blade. When the 
collocation is inside the box considered, the value 
of Sy must be calculated by the method of Ref. 18. 

For the curved wake 


n«4d 2 
- / 


f e 


v'(s'-s;) 
-* 7T^~ 


-ix' r' 
w 


B(- 


n j- d 2 


..3 


( 1 - ) + k , 2 6 2 z ’ 2 ] ds’dn' 

r* 
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-)[(l+i<’r;) 


(26) 


where r’ - [(x'-x ') 2 + B^y’-yl ) 2 + B 2 z’ 2 ) 1/2 , 
w m n 

x f * y 8 + y* sin 8 , y* m y j cos 8 , z 1 = d , 0 / 2 v, 

ds' » (dx , 2 +dy ,2 ) 1/2 - (y 2 +yj 2 + 2 yyj cos 6) 1/2 de, 

and d (~Rd f ) , the downward displacement of the wake 
per revolution, is assumed to be small. T^'s are 

evaluated numerically at the j ’ th spanwise strip by 
taking small Increments of 6 and n. 


Hovering Flight (Low Inflows) 


For hovering flight, y - 0, s* * yjQ, and 
ds f « yjd 8 . expression for , Eq. (25), can 

be simply modified by replacing y with zero and the 
wake integration for the j’th strip 


Results and Discussion 

A rectangular rotor blade of R/c * 10 was cho- 
sen and the. blade was assumed to extend from 0.1R 
to R. For the computation of the airload coefficie- 
nts, a grid of thirty six rectangular boxes consist- 
ing of six chordwise and six spanwise stations were 
used. The convergence of the results was tested by 
taking grid sizes of 6 x 8 (chordwise x spanwise) and 
8 x 6 . Rigid mode shapes for flapping and twisting 
oscillations are assumed so that 

c “ ry' + ax' 

w'(Xp.yp) * fiR[a’ (u+yM-ifx^) + y'if(u+y^)] 

2 z* 

- £2R[ (p+y M-ifXp)a' * if(y+y * ) (^) (-^) ] • (29) 

It should he noted that the above equation is valid 
only for the blade position at tfj * ir/2. For hover- 
ing flight, y - 0 in Eq. (29). 

The airload coefficients for Mach numbers 0 
0.8 for several values of ^ and wake spacing of two 

chords were obtained with reference to the blade's 
quater-chord axis. In Figs. (2) thru (5), selected 
airload coefficients for slow forward flight (y - 
0 . 1 ) are compared with the results obtained for 
hovering flight using both a helical wake model and 
two-dimensional strip theory. For this particular 
comparison, the reference velocity in Eqs. (23) and 
(24) was taken as the relative local velocity, U, 
and the tip Mach number was 0.8. From these plots, 
one can conclude that strip theory predicts substa- 
ntially larger values for the airload coefficients. 
One of the most Important observations one can make 
is that the curved wake changes the chordwise load 
distribution in such a way that the center of pre- 
ssure shifts forward of the quarter-chord axis 
position (see Fig. 4). 
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Figs. (6) and (7) compare the results by seve- 
ral mathematical models used for the hovering fli- 
ght case. The airload coefficients are referred 
to the tip velocity (OR) and this choice was made 
to indicate the trends of spanwise load distribut- 
ion. The Circular Wake model representation resu- 
lts in a substantial saving in computational time. 
For example, to obtain the airload coefficients for 
one set of geometric and flight conditions using 
6x6 grid on the blade, the Circular Wake model took 
only 1.5 minutes of computing time on IBM 360/65 
while the Helical Wake model took 2.5 minutes. 
Although the Circular Wake model seems to indicate 
the general trends of the airload coefficients, one 
should use the full helical wake to compute the 
airload coefficients accurately. 

Some typical results for hovering flight using 
the Circular Wake representation, compared with the 
results of two-dimensional strip theory, Ref. 7, 
are shown in Figs, (ft) thru (11). The results for 
the curved wake are in good agreement with the re- 
sults for strip theory for the inner blade sections; 
however, the agreement is poor towards the tip. 

Figs. (12) and (13) show the variation with axis 
position of M^, conveniently ref erred to as pitching 
moment damping airload coefficient, at spanwise 
stations of 0.475R and 0.925R, respectively. From 
these results, one can conclude that the agreement 
between the curved wake results and strip theory 
is good near the quarter-chord position for M - 0 
and M - 0.8 but it becomes very poor as the axis 
is moved towards the trailing edge. 
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Fig. 1 Schematic Diagram of Rotor Blade and 
its Wake. 
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